Phenomenological theory of the giant magnetoimpedance of composite wires 
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Composite wires with a three-layered structure are known to show a particularly large magne- 
toimpedance effect. The wires consist of a highly conductive core, an insulating layer and an outer 
ferromagnetic shell. In order to understand the origin of the effect a theory based on a coupling of 
the Maxwell equations to the Landau-Lifschitz-Gilbert equation is suggested. The theory is phe- 
nomenological in the sense that it does not account for a domain structure. However, theoretical 
results nicely reproduce those obtained in various measurements. Furthermore, an upper limit of 
the magnetoimpedance ratio for a given combination of materials can be determined. 



The giant magnetoimpedance (GMI) effect refers to a 
huge change in the impedance upon the application of 
an external static magnetic field H to a ferromagnetic 
conductor. This phenomenon was already mentioned in 
the 1930s, when the behavior of nickel-iron alloy wires 
was investigatecUi 2 -. In the early 1990s the GMI effect was 
analyzed more thoroughly for cobalt-based amorphous 
wires with micrometer dimensions by several groups^—. 
The relative change of the impedance Z is measured in 
terms of the quantity AZ/Z, called the GMI ratio, which 
is defined as 



AZ = \Z(H )\ 
Z 



\Z(H ref )\ 



Z(H re f) 



(1) 



where H re f represents a reference value of the applied 
magnetic field. The latter is usually taken to be its max- 
imum value at which the samples can be regarded as 
magnetically saturated. The ratio can easily reach up to 
several hundred percent and at the same time exhibit 
a strong field dependence in the low field regime. This 
makes the GMI effect promising for applications in sensi- 
tive magnetic field detection 7 ^. In order to increase the 
magnitude of the GMI ratio a wide variety of materials, 
material combinations and sample geometries have been 
analyzed 9 . - — . 

Heterogeneous structures, such as layered wires, are of 
particular interest^. Among them are composite wires 
consisting of a highly conductive core enclosed by a fer- 
romagnetic shel l 15 ! 16 . It has been found that the GMI 
ratio can be further enhanced by inserting an insulating 
layer between core and shell 1 — - — . Ideally, the presence of 
such an insulating layer restricts the flow of the driving 
AC current to the inner core only. 

In the present work we focus on establishing a model 
for the calculation and prediction of the GMI ratio in 
composite wires. The composite wires under investiga- 
tion are composed of a highly conductive nonmagnetic 
core, an insulating layer and a ferromagnetic shell. In 
order to model the impedance of the composite wires 
the distributions of the electric and magnetic fields are 
determined in each layer and an expression for the per- 
meability of the ferromagnetic shell is derived. 



The conductive core is characterized by its radius r c 
and its axis is along the z axis. When a driving cur- 
rent Ie~ lult with frequency / = uj/2-k flows through the 
core, the field distribution is found by solving Maxwell's 
equations in cylindrical coordinates neglecting displace- 
ment currents^ and by utilization of Ampere's law as a 
boundary condition. This restricts the driving current to 
the core and accounts for a perfect insulation of the inner 
core from the ferromagnetic shell by the interlayer. The 
amplitudes E z of the axial electric field and H v of the cir- 
cumferential magnetic field are represented by functions 
that solely depend on the radial coordinate r and can in 
CGS units be expressed as 20 



E (r) = C x J °( fccr ) 
J\{k c r c ) 



H v (r) = H T 



Ji(fc c r) 
Ji{k c r c ) 



(2) 
(3) 



with Hi = 2I/cr c and C\ = Ik c /2i:a c r c . Here a c 
denotes the conductivity of the core material and c the 
vacuum speed of light. J n represents Bessel functions of 
the first kind and order n. k c is defined by k c = (l+i)/S c , 
with the skin depth S c = c/ \/2iruja c of the core. 

In the insulating layer of thickness t\ there is no elec- 
tric current present. The solution of Maxwell's equa- 
tions for the field amplitudes, which are continuous at 
the conductor-insulator interface, is then given by 



E z {r) = C 2 -C l k c r c \n(pj , 



2J 

HM = - , 
cr 



(4) 
(5) 



with C 2 = CiJ (k c r c )/ Ji(k c r c ). 

The solutions within the ferromagnetic shell, the thick- 
ness of which is denoted by t 2 , are somewhat more in- 
volved. Due to B = fiH, with the permeability tensor 
ft, the equations that determine the field amplitudes H v 
and H z become now coupled differential equations. As 
a consequence the solutions of H v (r) and H z (r) in this 



2 



region are not Bessel functions. We assume them to be The angle is determined by the ratio of the magnitudes 
given by series^ that correspond to the expansions of of the external field and the anisotropy field strength, Ho 
Bessel functions of the first kind for the decoupled case: and Hk- By solving Eq. (9) one is able to obtain 



HM = C 3 -+J2 («n + c n In (1)) (1)" , (6) 

n 

^)=£(^nln(^))(£)\ (7) 

n 

with R = r c + t\ + <2- Inserting this approach into 
Maxwell's equations yields recurrence formulae for a n , 
b n , c n and d n and determines the constant C3. The com- 
ponents of the electric field can be derived from H v and 
H z by using Maxwell's equations for an ohmic conduc- 
tor with conductivity a m ignoring displacement currents: 
V x H = 47rcr m /c E. 

In Eqs. (5) and (6), and thereby in the four sets of 
coefficients, only three combinations of the permeability 
tensor's elements \x nm (n, m = r,ip, z) appear. They are 

defined by fj, v = fi w - (X^/fJ-rr, (J-z = fJ-zz - (J-rz/^rr 

and fx c = (J-<pz — HrzfJ>rip/ fJ>rr- The form of these elements 
of the permeability tensor have to be derived. This can 
be achieved by linearizing the Landau-Lifschitz-Gilbert 
equation 2 ^ 

d t M(t) = - 7 M(t) x H e// + -^-M(t) x d t M(t) . (8) 

Ms 

This equation describes the dynamics of the magneti- 
zation vector M(f) under the influence of the effective 
magnetic field H e ff- 7 denotes the gyromagnetic ratio 
of the electron, a the Gilbert damping parameter and 
Ms the saturation magnetization. Linearization can be 
achieved by subdividing the magnetization vector M(t) 
and the effective magnetic field H e ff into their equilib- 
rium parts M e<? and and into contributions Me" 1 "' 
and He~ lut that are excited by the driving current 23 . 
By using the relation M e? x H^L = and by omit- 
ting second-order terms in the excitations, the Landau- 
Lifschitz-Gilbert equation can be simplified to 24 i 25 

-M = M x frs% f + ^M eq ) + M eq x H . (9) 

This set of equations establishes a linear relation between 
the magnetization vector M and the magnetic field H 
through the susceptibility tensor \. The permeability 
tensor /} is then given by /t = 1 + 47r%. 

The equilibrium part H^jL is given by the sum of the 
external field Ho, which is applied along the z direction, 
and of a circumferential anisotropy field H^ 2 ®. Con- 
tributions originating from the exchange interaction and 
the demagnetizing field are not taken into account. The 
angle <p between the equilibrium magnetization and the 
wire axis is derived by minimization of the free energy. 
This results in 

H sin(0) - H K sin(0) cos(0) = . (10) 



fl v = 1 + H e ff COS 2 (» , (11) 

fx z = 1 + /i e// sin 2 ((/>) , (12) 
fi c = -fJ-eff sin(0) cos(0) . (13) 

Here, \i e j f denotes an effective permeability. For exter- 
nal field strengths sufficiently far above Hk (cos(</>) rj 1) 
the circumferential permeability fi v almost equals fi e ff- 
fi e f f consists of a real part Mg// an d an imaginary part 
p!' e tf and is defined by 

um(um + u K - lUOi) 
MS = 7 ; ; w : \ 2 ' ( 14 > 

with ujm = 'lir'yMs, ujk — jHk and ljq = r yHoCOs(4>)^- 
As a realistic example we assume 7 = 1.76 • 10 7 G _1 .s _1 , 
a = 0.1 and M s = 640 G. The field dependence of both 
the real and the imaginary part of the effective perme- 
ability is shown in Fig. [1] for two different values of Hk 
at a frequency of 2 MHz. The imaginary part exhibits 
a sharp peak at Ho = Hk, while the real part shows a 
double peak in this field regime. The permeability and 




H (Oe) 

FIG. 1. Dependence of the real and imaginary parts of the 
effective permeability on the applied magnetic field for Hk = 
and H K = 2 Oe at 2 MHz. 

field distribution previously derived are used to obtain 
the impedance Z of the composite wire by consideration 
of the total energy loss. The input power is given by 
P = ZI 2 . On the other hand, the energy dissipated per 
unit time from a volume element corresponds to the in- 
tegral of the Poynting vector S = c/47r E x H over the 
volume's surface. For a composite wire of length / this 
yields 

z = 4p £(- 1 ) nr »( E ( r *) x H M, ' ( 15 ) 

with r n defining the radial coordinates of the three lay- 
ers' surfaces. The field amplitudes of E and H are ob- 
tained from Eqs. (2) - (7). In the ferromagnetic shell the 
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field amplitudes explicitly contain combinations of the 
permeability tensor's elements with their strong field de- 
pendence. Ultimately, the impedance Z is a function of 
the amplitude I and frequency / of the driving current, 
of the wire geometry and additionally, through the per- 
meability tensor, of the applied magnetic field and the 
anisotropy field strength. 

With this expression for the composite wire impedance 
Z the dependence of the GMI ratio on the external field 
strength and the frequency of the driving current can be 
modeled. The GMI ratio exhibits a distinct maximum 
close to the magnitude of the anisotropy field. This fact is 
also observed in many experiments^. Beyond this value 
Hk, the GMI ratio starts to decrease. This yields that in 
case of a ferromagnetic layer not showing any magnetic 
anisotropy one would expect a maximum of the GMI ra- 
tio at Hq = 0. 

In the calculations we assume a typical experimen- 
tal scenario: r c — 30 /im, t\ = t% = 0.1 r c . The re- 
maining parameters are taken to be <r c = 5 • 10 17 s — 1 , 



10 



16 



H 



ref 



80 Oe and I = 10 mi. Here 



cr m represents the conductivity of the ferromagnetic shell 
material. Figure [5] depicts the modeled curves of the field 
dependence of the GMI ratio for two different values of 
Hk and two different frequencies of the driving current. 
It is obvious that the frequency of the driving current 
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FIG. 2. Field dependence of the GMI ratio at 200 kHz and 
2 MHz for Hk = and Hk = 2 Oe modeled for a composite 
wire of a 30 fim conductor core, a 3 /im insulating layer and 
a 3 //m ferromagnetic shell. 

has a large influence on the magnitude of the GMI ratio. 
The frequency dependence of the GMI ratio at different 
thicknesses of the insulating layer is shown in Fig. [3] 
The maximum of the GMI ratio at Hq ps Hk strongly 
depends on the thickness of the insulating layer with re- 
spect to its magnitude but not with respect to the fre- 
quency at which it occurs. The present theory predicts 
maxima in the 100 kHz range. This fact is corroborated 
by experimental results obtained for samples which are 



comparable to those assumed her o 18 ' 28 . Our results also 
predict an upper limit of the GMI ratio for t± — > and 
all other parameters as mentioned before. For a vanish- 
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FIG. 3. Frequency dependence of the GMI ratio of composite 
wires with a 30 /im core, a 3 fim ferromagnetic shell and 
different thicknesses of the insulating layer. 



ing insulator thickness AZ/Z approaches 1000%. In this 
case the vanishing insulator layer thickness is assumed 
to preserve the restriction of the driving current to the 
core region. This limit expresses that predicted GMI ra- 
tios are the highest for insulator layers prepared as thin 
as possible. If furthermore ti — > oo is assumed, even 
3000% are reached. These upper limits, which can easily 
be calculated for any materials combination, provide an 
important guideline for experimental research. 

In conclusion, a model to calculate and predict the 
GMI ratio of composite wires in a very general form has 
been established. The impedance of composite wires with 
insulating interlayer was calculated by separately solv- 
ing Maxwell's equations in each layer. The solution in 
the ferromagnetic layer additionally required the deriva- 
tion of an expression for the permeability tensor. This 
was achieved by linearization of the Landau-Lifschitz- 
Gilbert equation which involves several approximations. 
The field dependence of the GMI ratio for different values 
of Hk and the frequency dependence for different insu- 
lating layer thicknesses have been modeled. For a given 
frequency the GMI ratio reaches its maximum at an ex- 
ternal field strength equal to the value of the anisotropy 
field. The frequency spectrum indicates large GMI ra- 
tios in the 100 kHz range. The model predicts high GMI 
ratios for thin insulating and sufficiently thick ferromag- 
netic layers. The obtained values of 1000% for very thin 
insulating layer thicknesses can be compared to exper- 
imental results, where maximum GMI ratios of 800% 
have been reported^. It is interesting to mention that 
for other materials theoretical GMI ratio maxima above 
3000% were predicted 2 ^, which is in the same order of 
magnitude as our results. 
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